Discrete and surface solitons in photonic graphene nanoribbons 
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We analyze localization of light in honeycomb photonic lattices restricted in 
one dimension which can be regarded as an optical analog of ( "armchair" and 
"zigzag") graphene nanoribbons. We find the conditions for the existence of 
spatially localized states and discuss the effect of lattice topology on the proper- 
ties of discrete solitons excited inside the lattice and at its edges. In particular, 
we discover a novel type of soliton bistability, the so-called geometry-induced 
bist ability, in the lattices of a finite extent. 
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The studies of a monolayer of graphite sheet, called 
graphene, have attracted growing attention due to many 
interesting transport properties of electrons^. More- 
over, semi-infinite graphene and finite stripes of graphene 
(called graphene nanoribbons) with zigzag edges support 
peculiar electronic states with nearly flat dispersion. 

The interesting phenomena in graphene structures are 
not limited to the electronic systems, and they have di- 
rect analogs in the physics of photonic crystals^} and 
photonic lattices^. As a matter of fact, many of the 
phenomena are generic to honeycomb lattices and can ap- 
ply to electromagnetic waves in photonic lattices, quasi- 
particles in graphene, and cold atoms in optical lattices. 

All the problems considered for electronic properties 
of graphene are linear, and no nonlinear effects were dis- 
cussed so far. However, the photonic analogy suggests 
not only the study of nonlinear effects in graphene-like 
structures such as spatially localized nonlinear mode*^, 
but also a possibility of direct experimental verifications 
of many of the predicted phenomena, f or bo th hexagon 
and honeycomb two-dimensional lattices^S 

In this Letter we employ the analogy with graphene 
nanoribbons, and study localization of light in honey- 
comb photonic lattices of a finite extent, an optical ana- 
log of graphene nanoribbons. We find the conditions for 
the existence of spatially localized states and reveal the 
substantial influence of the lattice topology (i.e. "arm- 
chair" or "zigzag" ) on the properties of discrete solitons 
excited inside the lattice or at its edges. 

We consider a two-dimensional honeycomb photonic 
lattice with a finite extension in one dimension, an op- 
tical analog of the graphene nanoribbons. Such pho- 
tonic stripes can have two distinct geometries, which can 
be classified by employing the graphene terminology as 
"armchair" and "zigzag" structures, as shown in Figs. 1 
(a,b), respectively. In the framework of the coupled- 
mode theory, the electric field £ (r) propagating along 
the waveguides can be presented as a superposition of 
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Fig. 1. Schematic of optical lattice nanoribbons with (a) 
armchair and (b) zigzag graphene geometries. 



the waveguide modes, £(r) = ^n^n<K r — n ), where S n 
is the amplitude of the (single) guide mode 0(r) centered 
on site with the lattice number n = (711,712)- The evo- 
lution equations for the modal amplitudes E n take the 
well-known form, 



d£. 



dz 



+ V J2 £m+7|£n| 2 £n=0 



(1) 



n 1 ,n 2 



where n denotes the position of a guide center, and V is 
the coupling in the lattice. The nonlinear parameter 7 is 
normalized to 1 for the focussing nonlinearity. 

Next, we analyze the stationary localized modes of 
Eq. (fTl in the form £ n (z) = i? n exp(z/3z), where the am- 
plitudes E n satisfy the nonlinear difference equations, 



- pE n + V J2 E ™ + 7|£n| 2 £n = 



(2) 



n 1 ,n 2 



We consider a nonlinear case for which the linear 
regime can be achieved in the limit P — >• where P = 
^2 n \E n \ 2 is the mode power. For a given value of /?, the 
system of stationary equations pi is solved numerically 



by a multidimensional Newton- Raphson scheme. As we 
are interested in the modes localized inside the structure, 
we look for the localized solutions with the maxima near 
the center decaying quickly along and across the stripe. 
In order to visualize the field in the lattice, we present 
the field as U(x, y) = Y, n ,m C n,m(t>(x -n,y - m), where 
(j) is a guided mode of a single waveguide centered at the 
site (n, m). For the latter function, we assume a generic 
form, (j)(x,y) = exp[— (x 2 + y 2 )/cr], taking a = 0.1. Fig- 
ure ^[a) shows an example of the modes localized in the 
optical graphene stripe. 
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Fig. 2. Optical localized modes for the graphene nanorib- 
bon with the armchair geometry, (a) Three-dimensional 
intensity profile of a typical localized mode in a width 
3 nanoribbon (see Fig.2(c)). (b-d) Power vs. propaga- 
tion constant for the fundamental localized modes in a 
nanoribbon with the cross-width of (b) one, (c) two, and 
(d) four elementary cells. The insets show the corre- 
sponding structures with the shading of the intensity dis- 
tribution corresponding to the nonlinear mode at /3 ~ 3. 

We find that the results vary substantially depend- 
ing on the width of the stripe. In both cases, we find 
that relatively narrow stripes show the properties of one- 
dimensional nonlinear chains, where the spatially local- 
ized modes do not exist in the linear limit but split off 
the edge of the continuous spectrum for the power depen- 
dence P(/3). For P — »> 0, this curve approaches the value 
/? m that coincides with the edge of the linear band. In 
particular, for the narrow stripe of Fig. 2(b), we find an- 
alytically that the linear dispersion is described by two 
branches, /3i j2 (fc) = V[3 ± 2(1 + cos/c) 1 / 2 ] 1 / 2 , so that 
f3 m = V(3 + 2V/2) 172 which for V = 1 gives (3 m « 2.414, 
corresponding to the cutoff value in Fig. 2(b). 

For wider stripes, we observe the appearance of a kink 
in the power dependence and the corresponding mode 



bistability [see, e.g., Figs. |2jd)]. This kink will disap- 
pear for much broader (width 8) stripes, so the lattice of 
an intermediate extent demonstrates a crossover between 
one- and two-dimensional lattices. More importantly, the 
bistable dependence shown by the function P(/3) demon- 
strates the first, to the best of our knowledge, example of 
a geometry-driven bistability of solitons. Figure 4 shows 
stable modes on both sides of the bistable curve of Fig. 
2(d), obtained by dynamical evolution of the left mode 
of the branch (Fig. 4(a)), with j3 = 3.0, perturbed by 
adding an amount of power exceeding the maximum sup- 
ported by that branch. After some evolution time, the 
system falls into a mode of the right branch (Fig. 4(b)), 
characterized by a sharper localization, in view of its 
closer proximity to the anticontinuum limit. 

Surprisingly, the localized modes in the lattice with the 
zigzag geometry demonstrate a very different behavior 
with almost no crossover regime. Figures pfa,b) show 
the power dependencies for two types of "zigzag" stripes 
created of a honeycomb photonic lattice of a finite extent. 
In the weakly nonlinear regime the localized modes of 
narrow stripes do show the properties of one-dimensional 
discrete solitons similar to the modes in the armchair 
geometry. In particular, for the stripe of Fig. 3(a) the 
dispersion can be found in the form j3{k) = (V r /2)[1 + (1-|- 
16 cos 2 /c) 1 / 2 ], so that the cutoff value f3 m = (V/2)(l + 
yYJ) which for V = 1 gives /3 m ~ 2.56. For broader 
stripes we do not observe the crossover regime, and the 
power dependence acquire the genuine two-dimensional 
characteristics, see Fig. |3|b). 
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Fig. 3. Optical localized modes for the graphene nanorib- 
bons with the zigzag geometry, (a) Power vs. propaga- 
tion constant for the fundamental localized modes in a 
nanoribbon with the width of (a) one and (b) four ele- 
mentary cells. The insets show the corresponding struc- 
tures with the shading of the intensity distribution cor- 
responding to the nonlinear mode at /3 ~ 3. 

Finally, we analyze surface modes in such photonic 
structures. Existence of novel types of discrete sur- 
face solitons localized in the corners or at the edges 
of two-dimensional photonic lattices^^ have been re- 
cently confirmed by the experimental observation of two- 
dimensional surface solitons in optically-induced pho- 
tonic lattices^ and two-dimensional waveguide arrays 
laser-written in fused silic al 15 * 1 ^. These two-dimensional 
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Fig. 4. Dynamical transition between the two stable lo- 
calized modes on both sides of the bistable branches of 
Fig. 2(d). After vesting enough power on initial mode 
with small propagation constant (a), the system evolves 
dynamically to mode with larger propagation constant 

(b). 



nonlinear surface modes demonstrate novel features in 
comparison with their counterparts in truncated one- 
dimensional waveguide arrays^"^. In particular, in a 
sharp contrast to one-dimensional discrete surface soli- 
tons, the mode threshold is lower at the surface than in 
a bulk making the mode excitation easier^. 

Here, we employ our photonic nanoribbons and study 




Fig. 5. Examples of surface modes in a honeycomb opti- 
cal lattices. (a,b) Power vs propagation constant for for 
a nanoribbon of width 3 in the armchair and zigzag ge- 
ometries, respectively. The upper (lower) curve refers to 
a mode centered at a boundary site with two (three) near- 
est neighbors, as shown in the insets. The solid (dashed) 
curves denotes stable (unstable) portions while the ver- 
tical dashed line denotes the position of the linear band. 

localization of light at the edge of the lattice. We reveal 
that the effectively one-dimensional nature of the waveg- 
uide created in a two-dimensional lattice leads to the lo- 
calized surface modes which resemble the properties of 
the discrete surface solitons in a waveguide arrays^SELlD. 

Figures 5|a,b) show several examples of low-order non- 
linear surface modes, for both armchair and zigzag ge- 
ometries of the honeycomb photonic lattice, respectively. 
These modes do not have their linear counterpart and re- 
quire a threshold power for their excitation. The stability 
analysis of those surface nonlinear modes show that the 
well-known Vakhitov-Kolokolov stability criterium seems 
to hold, so that the branches with the positive slope in 
Figs. p[a,b) describe stable nonlinear surface states. 



In conclusion, we have studied localization of light 
in two-dimensional finite-size honeycomb photonic lat- 
tices, the so-called photonic graphene nanoribbons. We 
have revealed an important effect of the lattice geom- 
etry on the existence and properties of spatially local- 
ized modes and discrete solitons. We have demonstrated 
that the discrete solitons reveal an interesting feature of 
the geometry-induced bistability in the lattice of a finite 
width. Our results are generic to honeycomb lattices of 
a different nature, and they can apply not only to elec- 
tromagnetic waves in photonic lattices, but also to quasi- 
particles in graphene and cold atoms in optical lattices. 
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FB0824/2008), and by the Australian Research Council. 
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